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Abstract. We introduce a family of rotationally invariant random matrix ensembles 
characterized by a parameter A. While A = 1 corresponds to well-known critical 
ensembles, we show that X 1 describes "Levy like" ensembles, characterized by 
power law eigenvalue densities. For A > 1 the density is bounded, as in Gaussian 
ensembles, but A < 1 describes ensembles characterized by densities with long tails. 
In particular, the model allows us to evaluate, in terms of a novel family of orthogonal 
polynomials, the eigenvalue correlations for Levy like ensembles. These correlations 
differ qualitatively from those in cither the Gaussian or the critical ensembles. 
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1. Introdcution 

Gaussian Random Matrix Ensembles (RMEs) were proposed by Wigner about half a 
century ago to describe the statistical properties of the eigenvalues and eigenfunctions 
of complex many-body quantum systems in which the Hamiltonians are considered 
only in a probabilistic way P, [2]. Over the past several decades, they proved to be a 
very useful tool in the studies of equilibrium and transport properties of disordered 
quantum systems, classically chaotic systems with a few degrees of freedom, two- 
dimensional gravity, conformal field theory and chiral phase transition in quantum 
chromodynamics [3]. This wide applicability results from certain universal properties of 
the correlation of the eigenvalues known as the Wigner distributions, as opposed to the 
Poisson distributions that result from random sets of eigenvalues. 

More recently, attempts have been made to construct generalized ensembles that 
show a crossover from a Wigner to a Poisson distribution as a function of a parameter, 
as seen in many physical systems [1]. One such generalization is the family of 'q- 
random matrix ensembles' (q-RMEs) [5l [6]. These were later shown to be models of 
'critical ensembles', with statistical properties different from the Gaussian RMEs and 
relevant for systems near a metal- insulator transition [7]. While there are other models 
that also describe critical statistics [8], one advantage of the q-RMEs is that they are 
rotationally invariant, and therefore can be analytically studied in great detail by the 
powerful method of orthogonal polynomials [2]. In particular, the differences between 
the Gaussian and the critical RMEs can be traced to the differences in the asymptotic 
properties of the classical vs. q-orthogonal polynomials [9l [10], illustrating how the 
universality of the Gaussian RMEs breaks down and gives rise to a different kind of 
universality for the critical ensembles. 

In this work, we take the generalization one step further to include 'Levy like' 
ensembles. Levy ensembles were introduced by Cizeau and Bouchaud (CB) [TT] 
where the matrix elements are drawn from a power law distribution according to 
P{Hij) ~ l/\Hij\^^'^, Hij S> 1, < /i < 1. The eigenvalue density for such ensembles 
falls off as l/x^"'"'^. Such matrices with a broad distribution of matrix elements have 
been studied in the context of a wide variety of systems including financial markets, 
earthquakes, scale free networks, communication systems etc. as well as quantum 
chaotic systems such as the Coulomb billiard and kicked rotor with singularities |12j . 
The statistical properties of the eigenvalues and eigenvectors of such matrices are, in 
general, quite different from the universal properties of the Gaussian or the critical 
ensembles. Indeed, it has been argued [H] that matrices of the Levy type are relevant 
for describing localization transition of interacting electrons in infinite dimensions. In 
particular, while the eigenvectors of the Gaussian RMEs are all extended, the Levy type 
matrices contain signatures of mobility edges, separating localized states from extended 
states within the spectrum. 

However, the CB model is very difficult to study because it is not rotationally 
invariant. While eigenvalue densities for certain Levy ensembles have been obtained 
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analytically [13] , little progress has been made in any systematic study of the eigenvalue 
correlations where novel universal features can be expected. In the present work we 
will consider 'Levy like' ensembles, characterized by a power law eigenvalue density 
l/x^^'^, without specifying whether such a density arises from a distribution of the 
matrix elements as in the CB model. For this purpose, we introduce a rotationally 
invariant model characterized by a parameter A. We show that A = 1 corresponds to 
a critical ensemble, while A > 1 and A < 1 describe Levy like ensembles with bounded 
and unbounded (long tail) densities, respectively. In particular the model allows us to 
calculate the two-level kernel for Levy like ensembles, from which all correlations can 
be evaluated, in terms of a novel family of orthogonal polynomials. These correlations 
turn out to be qualitatively different from either the Gaussian or the critical ensembles. 

Consider the set of all x Hermitian matrices M (we will restrict ourselves to 
the unitary ensembles only) randomly chosen with the following probability measure 
P^{M)dM (X ex^[—Tr{y{M))]dM, where V{x) is a suitably increasing function of x, 
Tr is the matrix trace and dM the Haar measure. By going over to the eigenvalues- 
eigenvectors representation, it can be shown that the joint probability distribution of 
the eigenvalues X = {xi, i = 1,2, . . . N) of the matrices can then be written in the form 

m 

N N 

P;(X)cx n ix.-x,rlle-'^^-'\ (1) 

l<i<j<N i=l 

Here the factor Y[{ 

from the Jacobian of a change of variables. 
The level correlations can be determined exactly by recognizing that the distribution 
can be written as a product of Vandermonde determinants of a set of (monic) 
polynomials (pni^) that are orthogonal with respect to the weight function w{x) = e~^^^^ 

oo 

j e-^(^Vn(a;)0m(a;)c/x = Smn- (2) 

— oo 

The main quantity of interest is the large N limit of the two-level kernel 

N-l 

Kl{x, y) ^ e-(^(^)+^(^))/^ Mx)My), (3) 

n=0 

from which all correlation functions can be obtained. The Gaussian RMEs follow when 
one chooses the 'confinement potential' V{x) = that defines the Hermite polynomials. 
Since all 'Freud-type' weight functions e~^^^\ with monotonically increasing polynomial 
V{x), lead to orthogonal polynomials with qualitatively similar asymptotic behavior in 
the large limit, they all share the same correlations as the Gaussian REMs. This is at 
the root of the wide applicability of the Gaussian RMEs. On the other hand, q-RMEs 
follow when one chooses V{x) ~ In^x for large x [5], which leads to the q-polynomials. 
All q-polynomials have similar large A^ asymptotic behavior, which differ qualitatively 
from those of the Freud like classical polynomials. Thus all q-RMEs characterized by 
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different q-polynomials share the same correlation properties as those of the critical 
ensembles, which are different from those of the Gaussian ensembles. 

2. Model for Levy-like ensembles 

It turns out that all q-RMEs have density of eigenvalues falling off as 1/x |14j . 
Comparing this with the fact that Levy-like ensembles have density falling off as l/x^~^^, 
this suggests that we consider a generalization of the q-RME model where the asymptotic 
In^ X behavior of the potential is extended to include other powers of the logarithm. We 
therefore introduce a one parameter generalization of the q-RMEs characterized by 

V{x) = —^[smh~'x]'+'; A > 0; q < 1. (4) 
ln(l/g) 

For A = 1, this is a model for critical ensembles [B], defined by the Ismail-Masson q- 
polynomials [10]. Since the asymptotic behavior of the potential is V{x) ~ ln^''"^(x), 
any A 7^ 1 is qualitatively different from a q-ensemble. However, there is no direct 
way to obtain orthogonal polynomials for any arbitrary non-trivial weight function. In 
particular, the orthogonal polynomials corresponding to A 7^ 1 are not known, and we 
can not write down the two-level kernel (Eq. 3) directly. We therefore follow an indirect 
but systematic method [15] that allows us to obtain the polynomials recursively for any 
given potential. It is well known that every orthogonal system of real valued polynomials 
satisfy a three term recursion relation [9], [10] 

X(f)n{x) = (l)n+l{x) + Sn4>n{x) + i?„0„_i(x). (5) 

Following Ref [15], we define a set of integrals 

/oo 
x'"e-2^(^Vn(x)rfx, (6) 
-CXD 

which in turn satisfy the recursion relation 

Qn,m Qn—l,m—l Sn—lQn—l,m Rn—lQn—2,m- (7) 

Thus, the determination of the coefficients i?„ and Sn necessary to calculate the 
polynomials of degree n < N — 1 requires only the knowledge of the 2N + 1 integrals 
go,m form = 0,1, ■••27V. 

3. Results 

We first note that for the present model, because we choose the potential to be 
symmetric, all Sn = 0. Therefore, Rn completely determines the polynomials. For 
simplicity, we will use a fixed value of the parameter ln(l/g) = 7 = 0.5 except in 
Figure [2] where we also show results for 7 = 1. Figure [1] shows the n-dependence of 
Rn for different values of A obtained by evaluating the Qo,m numerically. For all A, the 
large n behavior of Rn has the form 

RnOcq~-"''\ (8) 
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Figure 1. Log i?„ as a function of n for different values of A. Solid line corresponds 
to the critical ensemble A = 1, which separates two qualitatively distinct classes A > 1 
and A < 1. 



s 









7=0.5 ^ 
7=1.0 □ 




- a 












m 










m 

e 


s 

m 

t 


a 

a 















0.5 1 1.5 2 



Figure 2. The exponent a(A) as a function of A for two different values of 7 = ln{l/q). 
A = 1 corresponds to the q-polynomials describing the critical ensembles. 



In contrast, for all Freud like classical orthogonal polynomials, Rn oc n. Figure [2] shows 
the A dependence of the exponent a. 

For A = 1, a (A) = 1, which recovers the known recursion relatioijjl for the 
q-orthogonal polynomials. On the other hand for A > 1, a(A) < 1, while for 
A < 1, a (A) > 1. Note that the distribution of the zeros of the polynomials are 
determined by the exponent a, such that when a = 1, the logarithm of the zeros of the 
corresponding Ismail-Masson q-polynomials are uniformly distributed. For a < 1 the 
zeros are 'bunched together' while for a > 1 they are 'unbunched', highlighting their 
differences from the q-polynomials. Thus, A 7^ 1 defines a novel family of orthogonal 
polynomials that generalizes the q-polynomials. However, the asymptotic properties of 
these 'generalized q-polynomials' are not known. 

:|: Exact expression for all n is i?„ = (q^" — l)/4, see [TO] . 
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Figure 3. Density of eigenvalues for different values of A. A = 1 is the critical 
ensemble. For A > 1, the density has a sharp edge, while for A < 1 there is a true 
'Levy-like' long tail. 



The density of eigenvalues p(x) = K^{x, x) can now be obtained for different values 
of A from Eq. 3 by summing the products numerically. The results are shown in Figure O 
As already known, the density for A = 1 falls off as 1/x. For A 7^ 1, the density falls 
off as l/x^^^^'^\ as expected for Levy-like ensembles. However, for A > 1, /i(A) < 
and the density falls off slower than 1/x; the normalization condition (that there are 
TV eigenvalues in total for a matrix of size A^) then forces the density to have a sharp 
edge, similar to the edge of the semicircular density of the Gaussian RMEs. Thus, even 
though all A 7^ 1 have power law densities, the long tail characteristic of a Levy ensemble 
is cut off for A > 1 by a sharp edge. For A < 1, ;u(A) > and the density has true 
power law tails. It should be noted that nonextensive ensembles [I6] characterized by 
one parameter q also show a similar power-law behavior; for q > 1, the distributions of 
eigenvalue density show true long tails and for q < 1, the distributions have compact 
support. However, the parameter q in this case depends on the dimensionality N of 
the ensemble such that in the large N limit where universal behavior is expected, the 
maximum q allowed for the nonextensive ensembles approaches unity. 

Our results now allow us to obtain exact correlations for Levy like ensembles. 
Using Eq. 3 for the two-level kernel we calculate the 'unfolded' cluster functiorl§| 
Y{r) = |A']^(r)p where the variable r is such that the density is uniform and unity: 
p(r) = 1. It is only in this rescaled variable r that the universality of a given ensemble 
is revealed if it exists. The results are shown in Figure HI Note that cluster functions 
for A > 1 are intermediate between the Gaussian (V^) and the critical (^'") ensmbles 
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7 sm^vrrj 
_27rsinh[(7/2)r] 



(9) 



§ We restrict ourselves to a range of parameters where the kernel can be considered translationally 
invariant. 
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Figure 4. The cluster function Y{r) for Levy-like ensembles for different values of A. 
The dotted and the solid lines agree with the known results for the Gaussian and the 
critical ensemble (with 7 = 0.5), respectively. 



but for A < 1 it is qualitatively different, with peak positions of the lobes shifting towards 
smaller r values with decreasing A. This suggests that Y consists of a more general A- 
dependent argument as well as power of the sinh function that reduces to F*" in the 
limit A = 1. Note that our results are limited to finite system size N only, which makes 
any attempt to obtain the true asymptotic behavior of Y unreliable. Nevertheless, with 
sufficiently large N, we can clearly demonstrate that the correlation functions obtained 
from y for A 7^ 1 are qualitatively different from both the Gaussian and the critical 
ensembles. 

As an example of the correlation functions obtained from the cluster function, we 
evaluated the number variance within a range L shown in Figure [5l For A >> 1, 
the number variance tends to that of the Gaussian RMEs, as expected, due to the 
absence of long tails. On the other hand for A < 1, the number variance shifts towards 
the uncorrelated Poisson distribution. For A < 1, for large L seems linear in L, as it 
is for the critical as well as Poisson cases, although the slope depends on A in a nontrivial 
manner. We have also evaluated the gap function which, for critical ensembles, falls off 
slower compared to the Gaussian RMEs. We find that for A > 1 the gap function is 
in between the critical and the Gaussian RMEs, while for A < 1 it falls off slower than 
the critical ensemble. Again, finite system size prevents us from obtaining a reliable A 
dependence of the asymptotic behavior. Determination of the true asymptotic behavior 
of these correlation functions will require knowledge of the asymptotic properties of the 
novel orthogonal polynomials introduced here. 

4. Summary &: Conclusion 



In summary, we have introduced a rotationally invariant model of random matrix 
ensembles by defining a confinement potential V{x) that includes a parameter A. We 
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Figure 5. Number variance for Levy like ensembles for different values of A 

obtained the orthogonal polynomials associated with the weight function w{x) = e~^^^^ 
for different values of A by using a recursive method. For A = 1, we recover the 
well known q-polynomials that describe the critical ensemble, with eigenvalue density 
p{x) (X 1/x. For A 7^ 1, we obtain a novel family of 'generalized q-polynomials'. While 
p{x) oc l/x^^^*-^^ for all A 7^ 1, only A < 1 correspond to Levy like ensembles with true 
long tails. We use the generalized q-polynomials (for finite N) to evaluate the two-level 
kernel from which all eigenvalue correlations can be obtained. We find that for A 7^ 1, 
the correlations are of novel types, differing from both the Gaussian and the critical 
ensembles. Clearly, studies of the asymptotic properties of the generalized q-polynomials 
introduced above would be very useful in better understanding the properties of Levy 
like ensembles. 

Although we have considered only the variation with respect to the parameter A, 
the model has another parameter, q, which we kept fixed. We have confirmed that the 
two parameters are independent, so that a single effective parameter can not describe 
the separate dependence of the two parameters. It would be important to find out if 
the q-parameter can describe e.g. any "transition" within the Levy like ensembles, as 
argued in Ref. [11]. While the mobility edge or the localization transition has clear 
meaning in the context of electron transport in disordered systems, it would be very 
interesting if it exists in Levy like ensembles, in the context of complex systems like 
earthquakes and scale-free networks. 
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